% 给定参数
initial_height = 0; % 米
initial_velocity = 100; % 米/秒
gravity = 9.81; % 米/秒^2
time = 0:0.1:10; % 秒

% 计算没有空气阻力的情况
height_no_drag = initial_height + (initial_velocity * time) - (0.5 * gravity * time.^2);
velocity_no_drag = initial_velocity + (gravity * time);

% 计算有空气阻力的情况（阻力系数为 0.5）
drag_coefficient = 0.5;
height_with_drag = initial_height + (initial_velocity * time) - (0.5 * gravity * time.^2) - (0.5 * drag_coefficient * initial_velocity^2 * time.^3 / gravity);
velocity_with_drag = initial_velocity + (gravity * time) - (drag_coefficient * initial_velocity * time);

% 绘制曲线图
figure;
subplot(2, 1, 1);
plot(time, height_no_drag, 'b-', time, height_with_drag, 'r--');
xlabel('时间 (秒)');
ylabel('高度 (米)');
title('高度-时间曲线');
legend('没有空气阻力', '有空气阻力');

subplot(2, 1, 2);
plot(time, velocity_no_drag, 'b-', time, velocity_with_drag, 'r--');
xlabel('时间 (秒)');
ylabel('速度 (米/秒)');
title('速度-时间曲线');
legend('没有空气阻力', '有空气阻力');